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We investigate S=l antiferromagnetic quantum spin chain, whose exchange couphngs are 
strongly disordered. By the real-space renomalization group method, introduced by Ma, Das- 
gupta, and Hu, the renormalization group flows are analyzed numerically in a plain of the 
anisotoropy of the exchange coupling vs. the staggered magnetic field. As the result, the 
Heisenberg point, which has a zero average of the exchange coupling anisotropy, is specified as 
the unstable fixed point against the anisotropy. 
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§1. Introduction 

As is well known, in one-dimensional S=l antiferromagnetic Heisenberg (AFH) models, tiie 
ground state is quantumly disordered by the strong quantum fluctuation.^) There is a finite excita- 
tion gap above their non-degenerating ground state. Affleck et al. have exactly proved this fact for 
the model with biquadratic interactions, employing the valence-bond-solid (VBS) picture.^) The 
existence of the gap in the low-energy excitation makes the system robust against weak bond- 
randomness in the case with Jmin ^ where Jmin is the minimum bond strength in all and A is 
the magnitude of the excitation gap in the case without randomness. 

However, when the randomness becomes stronger, i.e. Jmin ^ the situation can be changed. 
Two far separated spins can form a singlet state even if they destroy the excitation gap of the 
bulk system. If the system contains an as large VBS cluster as the system size, the string order 
parameter still has a finite value. But if not, the string order parameter may vanish. This transition 
from quantum Griffiths (QG) phase to random singlet (RS) phase has been investigated numer- 
ically by the exact diagonalization method,^-* the density matrix renormalization group (DMRG) 
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calculation,^) the quantum Monte Carlo (QMC) simulation,^) and the real-space renormalization 
group (RSRG) method.^°) 

The RSRG method is convenient one to determine the explicit phase diagram in the ground 
state of a system with randomness. This method has been introduced by Ma, Dasgupta, and Hu to 
study the S=^ AFH chain,^) and developed to being applied to other systems both in the analytical 
way and the numerical oneJ~^^^ Among these works, we should point out a very important one 
by Fisher.^) He has fully studied the properties of the S=| antiferromagnetic chains with various 
types of random exchange coupling in the analytical RSRG method. Especially, he has concluded 
that the Heisenberg point is nothing but the unstable fixed point (XXX RS fixed point), at which 
there occurs a transition from the XX RS phase to the Z (Ising) antiferromagnetic (ZAF) phase in 
a random XXZ chain. 

As pointed out afterward, the RSRG method has the defect when it is applied to S>1 spin 
system. -"^^^ Therefore, we should take care to use it in such a way that few inadequate renormalization 
processes are involved. Some authors have created skillful procedures, such as the use of the 
effective Hamiltonian, 9-13, 15, 16) g^^^ keeping more degrees of freedom.^°) However, these skillful 
procedures cannot be applied in a fully consistent way to models with uniaxial anisotropy. The 
uniaxial anisotropy arises inevitably when the renormalization process is executed in the S>1 
quantum spin system with the magnetic field or the exchange anisotropy. Therefore, we should 
give up applying these skillful procedures to S=l random AF spin chain with the anisotropy. 

In the present work we therefore adopt the straightforward RSRG procedure and study numeri- 
cally the renomalization group (RG) flows around the Heisenberg point under staggered magnetic 
field. The latter very much reduces occurence of the inadequate renormalization processes, i.e., 
our present approach is effecient to examine the RG flows against the introduced anisotropy. It 
is found for the first time that the Heisenberg point, which has a zero average of the exchange 
coupling anisotropy, is specified as the unstable fixed point against the anisotropy. The result is 
in disagreement with the previous work by Saguia et al}^^ They studied almost the same model 
as we do, but without the staggered magnetic field, and concluded that the fixed point situates on 
the point with a finite average of the exchange anisotropy. A possible origin of the discrepancy is, 
we consider, the adequacy of the region where the decimation procedures have been executed. 

The construction of this paper is as follows. In §2, we introduce our model Hamiltonian. In §3, 
we breifiy show the RSRG method for our Hamiltonian and discuss the confidence of our numerical 
method. In §4, the numerical results are presented, and in §5 we give the summary of our analysis. 

§2. Model 

Our Hamiltonian for the strongly coupled pair of spins Si, S2 is given by 

Hq = JSi • S2 + LS1S2 
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-D^{Slf - D2{Slf - h^Sl - h2SI. (1) 
These spins are weakly coupled to the neighbors via 

n = KiSi ■ S'l + 1^282 • S'2 + MiSfS'i^ + M2S^S'i 

= H'- D^{S['f - D^{S'if - hsS[' - h^S'i. (2) 

Diagonalizing Hq, we obtain its cigen values, {Ek\k = 0,1,..., 8}, with their eigen states, 
{|/c) I = 0, 1, ... , 8}, where we take |0) as the two-spins ground state. As long as the magnetic 
field is staggered, i.e., /ii/i2 < 0, |0) is always non-degenerated. Introducing the magnetic field, we 
should have both the uniform anisotropy L and the uniaxial anistropy D even when we start the 
renormalization process from the state with all Ljj = and Di = 0. 

§3. Real-Space Renormalization Group Method 

In the present work we restrict our consideration to at T = 0. Taking into account H by the 
perturbation expansion and dropping 0{H^) terms, we can obtain the modified ground-state energy 
as 



Eo + {o\n\o) + ^ 



k=i ^o-Ek 



= E' + J'S[ ■ s'2 + L'S'i^S'i 

-D[{S['f - D'^iS'if - h[S[' - h'^S'i. (3) 

The all non-zero matrix elements of 

(0|7^|A;) = {Oin'lk) 

A; = 0,1,..., 8, (4) 

in eq. (3) should be required. After some algebra, we obtain the all recursion relations of these 

renormalized parameters as summerized in Appendix . 

Unfortunately, our recursion equations are so complicated that it would be difficult to be solved 

even in approximated forms. Therefore, we execute the decimation process numerically after our 

predecessors (see Fig. 1).^^"^^) Considering a simple case only with isotropic exchange couplings, 

we can easily find that this perturbative renormalization group for S>1 chain is not adequate.^^'^^^ 

In fact, we obtain the recursion relation for the isotropic exchange coupling of a spin-S chain as 

j,^ 2S{S + l)K,K2 
3 J 

within 0{Tt'^), and this coeffecient 2S(S-|-l)/3 is larger than unity for S>1. However, only in the 
strongly disordered limit with staggered field, the singlet-triplet energy-level differences are typically 
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so larger than the nearest neighbor couphngs that we can rarely find the inadequate decimation 
processes, in which J' > J. Actually in our numerical analysis such inadequate decimation processes 
occur within 3% to the whole decimation processes. We should take into account that our numerical 
results contain the intrinsic errors within 3%, and point out that our numerical conclusion is obtained 
within this accuracy. 



Di hi J Da h2 



D-i hi 
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s; 
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D'l hi 



s; 



J' 



Da h2 



Fig. 1. Decimation process for our Hamiltonian. Closed and open squares belong to A- and B-sub- lattices respec- 
tively. 



In the present work we analyze systems in the strongly disordered limit, i.e., the initial distribution 
of exchange couplings is chosen as P(0 < Jjj <!) = !. We introduce the random staggered 

field with the distribution as P(0 < hf < /i™"^) = l/Zi'"'^^ and p{-h""'-' < hf < 0) = 
where h™"""^ is varied, and A and B mean sub-lattices' indexes. Decimating until the maximum of 
Jij becomes less then the cutoff energy Q, a parameter which specifies the degree of decimation, 
we obtain averages of the renormalized parameters corresponding to fi. In order to obtain the 
reliable results, we should check that Jij is the largest energy-scale among those involved in the 
problem such as \Lij\, \Di\, and in each renormalization step. If, for example, we start our 
decimation step from the point with {6) > 0, where 6 = Lij/ Jij and (• • •) denotes the average over 
the remaining bonds in each decimation step, already for Q 0.3 we are faced to the case with 
5^1 (see Fig. 3 below), where our RG rule cannot work well. For {6) < —0.2, on the other 
hand, it turns out that inadequate processes occur more than 3%. Therefore, we should restrict our 
quantitative analyses within only decimation steps for which 0, stays a few tenths of maximum of 
initial Jij. For this purpose, we need to prepare such a large system that the averaged values can 
be obtained within small errors even if Q stays in such a large value. Our system initially contains 
2^^ = 32768 spins, and the number of samples examined is 100. We call the isotropic point with 
both {6} = and (h) = the Heisenberg point, where |(^)| = + \{h^)\- 

§4. Results 

In Fig. 2, we show a typical RG flow in the \ {h)\ — {S) plain. As seen in the figure, the staggered 
field is an irrelevant perturbation on the line, where {S) = 0. We can also show that RG flows 
accumulate on the Heisenberg point along this line as far as < <C fi. 

Taking the initial value of {h) as strictly, we could only obtain uncertain results on the flows 
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Fig. 2. Renormalization group flow in the decimation process starting from {5} = 0.000, h"^"'" = 0.0040. 

near the Heisenberg point because the RG flows are intrinsically unstable. However, setting small 
but flnite values of {h) first, we can observe the systematic RG flows near the Heisenberg point as 
shown in Fig.3. In the figure, we show the RG fiows under decreasing J7 from the initial points with 
\{h)\ = 0.0004 and {S) ^ 0. The situations are completely different between (S) > and (6) < 0. 
The RG fiows approach toward the point (6) = —1 in the area (S) < 0, while they become the 
infinity as {S) — oo in the area (6) > 0. These two fixed points, (S) = — 1 and (S) — oo have 
been called the XX RS and ZAF point, respectively, and they correspond to the different phases. 
We can also confirm that the averaged single-ion anisotropy, (D), behaves differently between in 
the XX RS and in the ZAF phases. In the XX RS phase, (D) < 0, i.e., preferable to the low-spin 
state, and their absolute values are comparable to the cutoff energy J7. On the other hand, in the 
ZAF phase, (D) > 0, i.e., preferable to the high-spin state, and their absolute values evolve toward 
infinity as J7 being decreased. All these RG fiows branch off at the line (S) = 0. According to our 
more detailed analysis, no critical line can be found more than (S) = 0.000(3). 

Putting together our results discussed so far, we propose the schematic RG fiow diagram in the 
strongly disordered limit as shown in Fig. 4, in which one sees the Heisenberg point is the unstable 
fixed point. Any definite informations about the feature of the Heisenberg point cannot be obtained 
due to the intrinsic numerical difficulties in our formalism. However, referring to the results by 
other methods,^"^' ^'^^ we can identify this unstable fixed point as the XXX RS point. The feature of 
the Heisenberg point should be qualitatively the same as the Heisenberg point of the S=^ chain.^^ 
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Fig. 3. Renormalization group flow in the decimation process starting from h"^'^^ = 0.0004 and 6 = ±0.005, ±0.035. 



|<fA|+|<h^| 




Fig. 4. Schematic renomahzation flow diagram in the strongly disordered hmit. 



§5. Summary 

We have analyzed the strongly disordered S=l AF quantum spin chain by the numerical RSRG 
method. The RG flows around the Heisenberg point (S) = can be investigated by introducing the 
small staggered field, which is the irrelevant perturbation for our renormalization process. We have 
carefully executed numerical decimation processes in such a way that few inadequate processes are 
involved. Consequently, we have obtained the reliable RG flows around the Heisenberg point within 
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the sufficient accuracy. It is concluded that the Heisenberg point is the unstable fixed point against 
anisotoropy, and that it turns from(to) XX RS phase to(from) ZAF phase due to infinitesimal 
changes of anisotropic couplings. 
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Appendix: Recursion Relations 

The recursion relations for the renormalized parameters in eq. (3) are described as: 
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L' 



Eq 
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(A3)' + {B^f 
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E3 

Es — Eq 
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El — Eq 

\2 



- Eq 
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A2B2_ 

Eq 

^2 
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£Lil — rLiQ -C/2 ~ -C'O ^3 ~ ^0 -C'4 ~ -C'O -^'5 ~ ^0 ^6 ~ ^0 



h'l = h3-AQ 

h2 = hs — Bq 



(^3 



E3 — Eq 



+ 



+ 



iA4f 
-E4 — Eq 

{B,f 



{AqY 



E5 — Eq Eq — Eq 



{B5 



(Bf, 



Ez — Eq E4 — Eq £5 — Eq Eq — Eq 



(A-l) 



Hereafter, we use the abbreviations, L = L + Di + D2, h = hi — h2, hf = hi ± Di, and 
hf = h2 D2. In eqs.(A-l), Eq, Ei, and E2 should be represented differently, whether /i = or 
not. When ^ = 0, 



Eq 

El 
E2 



(J + L)2 + 8J2 + J + L 

J-L, 

^J{J + Ly + 8J^ - J -L 



(A-2) 



When h^O, 
2 



En 



/' 



pcos 



arccos 



P 



-3 



(J + L){(J + L)2 + 9j2-9/i2}] +27r(n + l)}') - J-L 



(J + L)2 + 6J2 + 3/i^ 



1/2 



n = 0,1,2. 



(A-3) 
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On the other hand, E^jE^, E^, and Eq can be described independently of h: 

Ei = ^{h^ + +x~), 

E5 = ^i-ht -ht -x+), 

Ee = ^i-ht-ht + x+), 

using the abbreviation x^ = [4J^ + {hf — /if )^]^/^. and also depends on h. 

When ^ = 0, 
Ao = Bo = A2 = B2 = 0, 
Ai = V2J{Ki + Mi)r^\ 

Bi = -V2J{K2 + M2)rQ\ 

As = [2J2 - {Eo + J + L){hl- + x-)]Ki{V2rors)-\ 
Bs = [2{Eo + J + L)-{hi-h2+ x-)]JK2{V2ror3)-\ 
^4 = [2J^ - (£^0 + J + Z)(/ir - - x-)]Ki{V2ror^)-\ 
Bi = [2{Eq + J + L)-{h-^-h2- x-)]JK2{V2rori)-'^, 
^5 = [2(£^o + J + L) + ht-ht- x+] Jifi(V2ror5)-\ 
^5 = [2J2 + (^0 + J + Z)(/i| -ht- x+)]K2(\/2ror5)-\ 
^6 = [2(^0 + J + L) + ht-h+ + x+] JKi(\/2ror6)-\ 
^6 = [2J2 + {Eq + J + L){ht -ht + x+)]K2{^/2rQr^)-\ 
ro = [{E^ + J + Lf + 2J^]y\ 

And when h ^ 0, 
Ao = -AhJ'^iEo + J + L){Ki+ Mi)ro-2, 
So = Ah.fi {Eo + J + L){K2 + M2)rQ^, 
Ai = -2hj\Eo + Ei + 2J + 2L){Ki + Mi){rori)-\ 
Bi = 2hfi{Eo + Ei + 2J + 2L){K2 + M2)(rori)-S 
A2 = -2hfi{Eo + E2 + 2J + 2L){Ki + Mi)(ror2)-\ 
B2 = 2hfi{Eo + E2 + 2J + 2L){K2 + M2)(ror2)-S 

As = {Eq + J + L- h)[2fi -{E0 + J + L + h){hi - /12 + x-)]Ki{V2ror3)- 
B3 = {Eo + J + L + h)[2{Eo + J + L -h) - {h^ ~ + x-)]JK2{V2rors)-^ 
A4 = {Eo + J + L- h)[2fi -{E0 + J + L + h){hl - /12" - x-)]Kx{V2rQri)- 



Bi = {EQ + J + L + h)[2{EQ + J + L-h)-{h^-h^- x-)]JK2{V2ron)~^, 
A5 = {Eo + J + L + h)[2{Eo + J + L-h) + ht-ht- x+]J Ki{V2ror^)-^ , 

B5 = iEo + J + L- h)[2J^ + {Eo + J + L + h){ht - - x'^)]K2{V2rQr^)-\ 

Aq = {Eq + J + L + h)[2{Eo + J + L-h) + ht -h+ + x+]JKi{V2ror6)-\ 

B6 = {Eo + J + L- h)[2.f + {Eo + J + L + h){ht-h+ + x+)]K2{V2rore)-\ 

r„ = [{{En + J + Lf- h^f + 2{{En + J + Lf+ J^] n = 0, 1, 2. (A-6) 

In eqs.(A-5) and (A-6), rs, r^, r^, and tq are common: 

r3 = [2x-{x- +h]: -h^)]^/^, 
r4 = [2x~{x'' - hi + /i2 )]^''^ 

r5 = [2x+(x+-/i+ + /i+)]i/2^ 

re = [2x+{x+ + ht -hi)]^/\ (A-7) 
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